Phase Bifurcations of Strongly Correlated Electron Gas 
at the Conditions of dHvA Effect 
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In a framework of catastrophe theory we investigate the equihbrium set for the system of strongly 
correlated electron gas at the conditions of dHvA effect and show that the discontinuities accom- 
panied the diamagnetic phase transition (DPT) is handled by Riemann-Hugoniot catastrophe. We 
show that applicability of the standard condition for observation of DPT a >1 where a is the 
differential magnetic susceptibility is valid only in the nearest vicinity of triple degenerate point 
corresponding to the center of dIfvA period, but for arbitrary value of magnetic field in every pe- 
riod of dHvA oscillations this condition is modified in accordance with the bifurcation set of cusp 
catastrophe. While at the center of dHvA period the symmetric supercritical pitchfork bifurcation 
gives rise to the second order phase transition on temperature, the deviation of magnetic field from 
the value corresponding to the center of dHvA period results in the change of the phase transition 
order from the second to the first one both on temperature and magnetic field. In the framework 
of developed theory we obtain good agreement with available experimental data. 
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I. INTRODUCTION 



The instability of an electron gas due to strong cor- 
relations at the conditions of dHvA effect resulting in 
diamagnetic phase transition (DPT) into inhomogeneous 
diamagnetic phase (IDP) with formation of Condon do- 
mains (CDs) [l], _2J is intensively studied both theoret- 
ically and experimentally 3-[IH- The realization of in- 
trinsic structure of IDP is governed by the competition 
between long-range dipole-dipole interaction and short- 
range interaction related to the positive interface energy 
with typical magnetic length of Larmor radius. The ob- 
served IDP in Ag by splitting of NMR line [2] was iden- 
tified as periodic domain structure with alternating in 
neighboring domains one-component magnetization. The 
further experiments on observation of IDP by method of 
/iSR spectroscopy revealed the arise of diamagnetic in- 
stability in Be, Sn, Al, Pb and In Q-Q, giving strong 
support to the idea of formation of CD structure [l| . 

The experimental demonstration of existence of IDP in 
normal metals 0,0- 01 is very convincing and support 
the basic ideas of the theory of DPT developed in [l3l |. 
but the full understanding of the properties of IDP is still 
lacking. The attempts of the quantitative analysis of the 
data on the basis of model of plane-parallel band patterns 
reveal the contradictions [8] , put questions on correctness 
of applicability of demagnetizing coefficient for adequate 
description of shape-dependent properties of IDP [J, P , 
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and raised important unanswered questions concern- 
ing the type of the diamagnetic ordering, irreversibility 
of DPT and morphology of the domain patterns. 

There are striking similarities between IDP in nor- 
mal metals and other strongly correlated systems which 
undergo phase transition on temperature and magnetic 
field with formation of complex macroscopic patterns, e. 
g. the type-I superconductors and thin magnetic films. 
The different technics, including the powder pattern and 
magneto-optic methods, successfully used for observation 
of intermediate state of type-I superconductors revealed 
a very rich structure which in spite of all it com- 
plexity, amazingly reminds the variety of domain struc- 
tures in thin film [l^. Recent experimental observation 
of CDs by a set of micro Hall probes at the surface of a 
plate-like sample of Ag demonstrates the formation of 
complex structure which differs from the expected reg- 
ular CD patterns separated by plane-parallel DWs [l| 
and calls for further development of the experimental 
techniques to study the IDP. The first experimental ob- 
servation of such an exotic phenomena as diamagnetic 
hysteresis by complex method including direct measure- 
ments by Hall probe, standard ac method with different 
modulation levels, frequencies and magnetic field ramp 
rates [4] allows to reconstruct the magnetization reversal 
in Be and confirms the possibility of the first order phase 
transition. The giant parametric amplification of non- 
linear response which proves to be sensitive method for 
investigation of phase transitions in superconductors was 
observed in a single crystal of Be in quantizing magnetic 
field , where the measured output signal amazingly fol- 
lows the shape of the DPT diagram and shows the asym- 
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metry and shift to the upper edge of the dHvA period in 
ascending magnetic field justifying the possibihty of the 
first order DPT and necessity of further development of 
the theory of DPT. 

Motivated by recent data on observation of DPT insta- 
bility 0-0], in the framework of catastrophe theory we 
investigate the equilibrium set for the system of strongly 
correlated electron gas at the conditions of dHvA effect. 
We show that in every period of dHvA oscillations the 
discontinuities of order parameter accompanied DPT is 
handled by Ricmann-Hugoniot catastrophe implying the 
standard scenario for the transition, e.g. DPT is of the 
second order at the center of dHvA period, weakly first 
order in the nearest vicinity of this point and is of the 
first order at the rest part of the dHvA period both on 
temperature and magnetic field. Thus, similar to other 
magnetic systems, e. g. the spin [l^ and metamagnetic 
ones flT'l, DPT can be realized in a number of different 
nonuniform phases depending on the shape of the sample. 

At r = K the fine, sub-band, structure of Landau 
level can results in DPT of the first order on magnetic 
field [3, but this effect based on treating the CD in- 
stability as an electron topological transition is negligi- 
bly small in comparing with the results obtained in the 
present publication. 

We show that the condition for DPT occurrence a =1 
where a = fio max{9Af/9i?} is a differential magnetic 
susceptibility is valid only in the nearest vicinity of the 
center of dHvA period, x —0, where x is an increment of 
magnetic field [13,]. This condition is violated for a; 7^ 
and is replaced by the generalized condition on a in accor- 
dance with the bifurcation set of cusp catastrophe. We 
discuss the conditions for realization of the DPT in the 
sample with taking into account the long-range dipole- 
dipole interaction and compare the results with available 
data. 
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FIG. 1: Reduced magnetization y as a function of inverse dif- 
ferential magnetic susceptibility 1/a at three different values 
of magnetic field x =0 (a), tt/36 (b) and 7r/3 (c). The solid 
(dash) lines correspond to the local minimum (maximum) of 
the free energy density G{y;a,x). Symmetric supercritical 
pitchfork bifurcation (a) of the system undergoing the second 
order phase transition at the center of dHvA period a; = 
for a = 1 (critical point) is replaced by the imperfect bifurca- 
tion diagram (b, c) with increase of | a; |< tt {2tt is a period 
of dHvA oscillations in reduced units) when the system is a 
subject for the first order phase transition. 



II. MODEL 

Equilibrium properties of strongly correlated electron 
gas at the conditions of dHvA effect in one-harmonic 
approximation can be described by Gibbs free energy 
G{y; a, x) = a cos {x + y) + , which is a functional of 
two conjugated variables x and y — —d^G [11]. Here, the 
small-scale field x = kfio{H — Ha) & [— tt, +Tr] is the incre- 
ment of the large-scale internal magnetic field fioH and 
applied magnetic field fioHa, y is oscillating part of re- 
duced magnetization, k = 2TrF/{^oHa)^ — 2tt/AH, F is 
the fundamental frequency of the dHvA oscillations cor- 
responding to the extreme cross-section of Fermi-surface, 
and AH is dHvA period. Minimization of the free energy 
G(y; a, x) with respect to y leads to well-known expres- 
sion for the magnetization [13] 

y — a sin 6 (1) 

with the magnetic flux density b = x + y being a func- 
tion of increment of internal magnetic field x and reduced 



magnetization y. Internal magnetic field x = xq + Xms 
depends on applied magnetic field xq and stray field Xms , 
which take into account the long-range dipole-dipole in- 
teraction dependent on the shape of a sample and has 
to be found by solving the Maxwell equations in mag- 
netostatic approximation by a self-consistent procedure. 
As it was suggested by Shoenberg [l^ , the magnetic flux 
dependence of oscillating part of magnetization y Eq. ([T]) 
at curtain conditions, when the differential susceptibil- 
ity a >1, results in arising instability of uniform phase 
(Shoenberg effect). At this condition both the magne- 
tization and magnetic induction become a multi-valued 
function of internal magnetic fleld x. Therefore, depend- 
ing on long-range dipole-dipole interaction at high mag- 
netic field and low temperature the system of strongly 
correlated electron gas can undergo the DPT with for- 
mation of IDP which corresponds to a global minimum 
of the free energy of the system. Due to the governing 
role of the long-range dipole-dipole interactions, the mor- 
phology of the DPT depends crucially on the shape of the 







FIG. 2: (color online) Shown are the equilibrium surface and 
bifurcation set (solid line) of cusp catastrophe for the system 
of electron gas under conditions of strong dHvA effect. At 
crossing the bifurcation set with change of control variables 
a and x, the system of strongly correlated electron gas is 
a subject for the DPT of the first order. The only point 
{a — 1, X — 0) where the system undergoes the second-order 
phase transition, e. g. the critical point, is a triple degenerate 
point shown by solid circle. 



sample which, in general, cannot be account for by the 
straightforward use of conception of demagnetizing coef- 
ficient [gj. 

The magnetization y Eq. ([1]) as a function of inverse 
reduced ampUtude of the dHvA oscillations 1 /a at three 
values of the increment of magnetic field x =0, 7r/3 and 
2tt/3 is plotted in Fig. [1] which illustrates the typical 
bifurcation behavior of the system, e. g. the existence 
of two stable states for magnetization at some values of 
the parameters a >1 and x, and the possibility of discon- 
tinues change of the order parameter through the disap- 
pearance of stable states when x ^ 0. 
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FIG. 3: (color online) Free energy density G — G{b) as a 
function of magnetic induction for fix value a =2.25 and three 
different value of x =0, 0.5 and 1.5 is shown for one dHvA 
period, (a) For x —0 and Oc =1, the system has twofold de- 
generate ground state: two minimums (circles) corresponding 
to the same value of G can give rise to the simple plane- 
parallel domain structure. (b,c) The increase of x followed by 
increase of breaks the equivalency of the steady states, re- 
sults in increase of the difference between the the energies of 
two equilibrium states and possibility of realization of variety 
of domain structures. For a =2.25 and x =1.5 the only min- 
imum implies the absolute stability of uniform diamagnetic 
phase. 



III. RESULTS AND DISCUSSIONS 
A. Condition of Diamagnetic Instability 

The discontinuities of the function b = b{a, x) where 
parameter a controls the amount of ordering, or the 
value of order parameter, and parameter x breaks the 
Z2 symmetry of order parameter, can be handled in the 
framework of catastrophe theory. The standard approach 
allows to obtain the canonical expression for the nor- 
mal form f{r];f3) — f3i — P2V + justifying the exis- 
tence of cusp catastrophe in the problem (Fig. Here, 
the new control variables j3i and (32 relate to the ori- 
gin ones (a, x) by means of relations (3i = —Q'^xa~^ and 
/32 =6~^(a — l)xa~'^ correspondingly, and rj —Qa^^b is 
the scaled phase-state variable. The direct calculations 
lead to the explicit form of the bifurcation set 



acos(Va^ - 1- I a; I) = 1. (2) 

defining the bifurcation curve on (a, x) plane of con- 
trol variables where the fold bifurcation occurs. Eq. ([2]) 
is a generalized condition for DPT occurrence. It can 
be simplified in the vicinity of triple degenerate point 
(a=l,x=0) resulting in standard cusp bifurcation set 

8(1 - af + = (-3-) 

which is a semicubic parabola. It follows from Eq. ^ 
that the usually accepted condition for DPT occurrence 
a =1 [l3| is justified only at the center of the dHvA 
period, x = 0, but violated for x 7^0 and replaced by 

a = flc , (4) 
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where Oc is defined by Eq. ([2]) is a field-dependent. In 
particular, for x — » 0, we obtain 
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.2/3 



(5) 



The character of bifurcation may be seen from the sur- 
face plot b = b{a, x), or the surface of stationary solutions 
of & = a; -I- a sin 6. The results of calculations are repre- 
sented in Fig. [2] which shows the equilibrium manifold 
of the system of strongly correlated electron gas for one 
period of dHvA oscillations in joint space of state and 
control variables (6, a, x) and the bifurcation set Eq. dJ) 
in (a, a;)-plane. At the center of dHvA period {x = 0) 
and additional condition a — > 1 the system exhibits the 
symmetric pitchfork bifurcation and undergoes the phase 
transition of the second order from homogeneous to inho- 
mogeneous state at decrease of temperature T (increase 
of a) with formation of CD structure the type of which 
is governed by the long-range dipole-dipole interaction 
and depends on the shape of the sample. The critical 
point a = l,a; = 0, e. g. triple degenerate point, is the 
only point where the DPT is of the second order. But, 
at a; 7^ 0, a discontinues jump of order parameter im- 
plies the possibility of existence of the first order DPT 
when the projection of equilibrium state on (a, a;)-plane 
crosses the bifurcation set with change of temperature T 
or magnetic field x (see, also Fig. [Ij . A simple graphi- 
cal illustration of above-mentioned effects is provided by 
corresponding potential function, G{b;a,x), drawn as a 
function of magnetic induction b for a = 2.25 and three 
different value of x in Fig. [3] which shows the possibil- 
ity of formation of CD structure due to existence of two 
minima around the center of dHvA period (Fig. [5^,,b) 
and the absence of the condition for realization of IDP 
for greater values of x when two minima coalesce and the 
uniform ground state is not degenerate (Fig.[31;). 




FIG. 4: (color online) (a) Magnetostatic field, Xma, and (b) 
internal magnetic fields, x, at crossing the bifurcation set from 
high field side are plotted as a functions of applied field, xo, 
for difi'erent values of demagnetizing coefficient, n. At given 
value of internal field, shown by horizontal dash line in (b), the 
values Xi {i =1-6) of applied field correspond to the increasing 
values of n. Both function, Xm3 and x, are odd functions of 
Xo, thus, only positive interval of xq is shown. 



B. Shape Effect 

The catastrophe theory allows to eliminate all the pos- 
sible stable equilibria (local minima) for a given pair of 
control variables (a, x) with the magnetic induction 6, as 
a phase-state variable, and establish the boundaries of 
CDs instability. Thus, the straightforward application of 
the results consists in evaluation of the range of existence 
of IDP, Axo , defined as a part of dHvA period occupied 
by CDs. The boundary between different phases is gov- 
erned by the bifurcation set of the cusp catastrophe. Due 
to the long-range character, the dipole-dipole interaction 
is sensitive to the shape of a sample which results in the 
dependence of the internal magnetic field x on experi- 
ment arrangement. The simplest way to account for the 
shape of the sample is the use of demagnetizing coefficient 
n for uniformly magnetized ellipsoid subjected to the uni- 
form magnetic field along principal axis. In this case, the 
magnetization and magnetic induction are single-valued 
functions of applied magnetic field, xq, the ground state 



of the system is uniform and the shape effects can be cal- 
culated numerically. This procedure allows to evaluate 
the conditions of occurrence of the diamagnetic instabil- 
ity. Of cause, the morphology of the domain patterns 
in the IDP cannot be account for by the use of the con- 
ception of demagnetizing coefficient, but has to be es- 
tablished by solving Maxwell equations. The results of 
calculation of the stray field, Xms, as well as the inter- 
nal field, X, as a function of applied field, xq, for several 
values of n when the uniform magnetization goes to the 
value defined by bifurcation set Eq. ^ are plotted in 
Fig. [5] which illustrates the well-known effect, e. g. the 
increase of n results in the growing difference between in- 
ternal and applied magnetic fields due to the increase of 
absolute value of stray fields. In particular, considering 
the samples with different demagnetizing coefficients, it 
follows from Fig.lDJb), that at other equal conditions, the 
phase boundary is reached at the higher values of applied 
field for the sample with the larger value of n. Therefore, 
IDP occupies the larger part of the dHvA period for the 
sample with n 1 which are preferable for observation 
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of CD instability. In connection to this, it is worth to be 
mentioned that the discovery of the Condon domains was 
made on a plate-Uke sample of Ag [2] and the diamag- 
netic instability was observed mostly by use the sample 
with the values of relevant demagnetizing coefficient in 
the range n « 0.5 - 1 (see, [H, jsj). 

The effect of "broadening" of IDP when the phase 
boundary is drawn in (Q;i,a;o) plane is illustrated by 
Fig. [5] which shows the first-order phase curves for several 
values of demagnetizing coefficient. In the calculation we 
neglected the overheating-undercooling effects. 

To compare the results of the theory with available ex- 
perimental data on observation of CDs instability we cal- 
culate the characteristics of IDP. The formation of IDP 
is characterized by two measured principal parameters, 
the average magnetic induction splitting, (56, defined as 
the difference, calculated at the center of dHvA period, 
between the values of magnetic induction in two adjacent 
domains and the range of existence of CDs, Sx, defined 
as a part of dHvA period occupied by non-uniform phase 
[lOj . Though the occurrence of IDP and the existence of 
CDs itself depends on long-range dipole-dipole interac- 
tion, the average value of magnetic induction splitting, 
6b, depends on the differential magnetic susceptibility 
and with reasonable accuracy is defined in a explicit form 

by El 
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6b — 2a sin — 0. 



(6) 



As far as concerned the range of existence of IDP 6x, it is 
a shape-dependent quality due to dependence of internal 
field on dipole-dipole interaction sensitive to the sample 
shape, as illustrated in Fig. Inasmuch as both char- 
acteristics of IDP originate from the same Eq. ([1]), they 
can be related to each other by eliminating a. It results in 
rather complicated functional dependence which is han- 
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FIG. 5: (color online) The phase boundary is shown for dif- 
ferent values of demagnetizing coefficient n in one period of 
dHvA oscillations. At given value of n, the inner (outer) re- 
gion, designated by I {II), corresponds to non-uniform (uni- 
form) diamagnetic phase. For fix a (horizontal dash line) the 
part of dHvA period occupied by IDP is wider for the sample 
with larger demagnetizing coefficient. 
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FIG. 6: (color online) Family of curves 56 — 5b{Ax, n) is 
shown for one dHvA period. The slope of the curves at Sx ^ 
is defined by inverse demagnetizing coefficient, 1/n, in accor- 
dance with Eq. (O justifying the Shoenberg theory of shape 
effects at the vicinity of critical point (a=l, a;=0). The sym- 
bols accompanied by reference number represent the experi- 
mental data. 



died numerically, but in a limit 6x, 6b 
form of linear function 

6b w —6x. 

n 



has a simple 
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From mathematical point of view Eq. ([7]) is fulfilled in 
an infinitely small vicinity of the triple degenerate point 
(x =0, a = 1), justifying the applicability of Shoenberg 
theory of shape effect at this critical point. Accounting 
for the shape of the sample we calculate the family of the 
curves 6b = 6b{6x, n) plotted in Fig. [6l The graphical 
representation of function 6b = 6b(6x, n) which relates 
the characteristics of IDP to the shape of the sample is 
convenient to use for comparing the theoretical results 
with experiment. For this, we introduce also the char- 
acteristics of IDP which can be measured (or evaluated) 
directly in experiment, 6H — 6x/k and 6B = 6x/k. 

The formation of IDP is accompanied by irreversible 
behavior of the magnetization curve which can be used 
for experimental investigation of IDP. The experimen- 
tal detection of diamagnetic hysteresis was performed on 
Be single crystal rod- like sample of size 8x2x1 mm^ ^ . 
The prolate ellipsoid which approximates the sample has 
demagnetizing coefficient n = 0.04. The hysteresis was 
observed at wO.25 part of the dHvA period which is 27r in 
reduced units. It is reasonable to associate the range of 
existence of IDP with the range of existence of hysteresis. 
Thus, we obtain 6x Ril.57. Unfortunately, the magnetic 
induction splitting is not reported in experiment, but it 
can be calculated by use of the value a due to Eq. 
At the conditions of experiment, T =1.3 K, Td —2 K 
(Dingle temperature) and ^oHa =3.6 T, the differential 
magnetic susceptibility calculated in the model of slightly 
corrugated Fermi surface ^ll| is a Ri2 which with the aid 
of Eq. ^ gives 5b w3.8. To verify the consistency of this 
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value of 5b we calculate the absolute value of induction 
splitting SB «78.6 G using the reported value of dHvA 
period AH =130 G and compare it with the experimen- 
tally detected global change of induction AB Ril21 G at 
one dHvA period. The inequality SB < AB confirms 
the reliability of calculated value of 6b «3.8. The point 
corresponding to the values of 6x «1.57 and 6b «3.8 is 
plotted in Fig. [5] and lies at the curve for n «0.04, as it 
is expected. 

The series of experiments on observation of CDs by 
fiSR spectroscopy p-Q have proved the occurrence of 
diamagnetic instability of strongly correlated electron 
gas in a number of normal metals, revealed the pos- 
sibility of the first-order phase transition by detecting 
the overheating-undercooling effects and actually served 
as one of the stimulation factors for development of the 
theory. As an example, we consider the data for Be [GJ 
which provide the value of magnetic induction splitting 
and range of existence of CDs with high accuracy of « ±1 
G in applied field of several Tesla. The /xSR measure- 
ments in were performed on plate-like samples of size 
1x1x8 mm^. At the conditions of experiment, T =0.8 
K and noHa =2.75 T, the experimental value of induc- 
tion splitting 6B = 28.8±1.4 G which gives 6b = 2.31 
with use of the reported dHvA period AH —78.2 G. 
The value of the range of existence of IDP is 6H «27 
G or 6x «2.2. The plotted in Fig. [H point, 6x «2.2 and 
6b «2.31, belongs to the curve with n «0.8 which is close 
to the theoretical value of n =0.75 for oblate ellipsoid 
approximating the sample. 

The famous experiment on observation of CDs instabil- 
ity by NMR technique on the plate-like sample of Ag 
of size 8x8x0.8 mm^ at helium temperature and applied 
field fioHa =9 T provides both values 6B —6H « 11 
which allows to calculate the reduced values 6b —6x « 
4.32 with use of experimental value of dHvA period 
AH =15.9 G. The corresponding point plotted in Fig.lH] 
lies at the curve parameterized by n w 0.78 which is close 
to the value n w 0.84 obtained for inscribed oblate ellip- 
soid approximated the sample. Thus, there is a good 
agreement of the theoretical results with the experimen- 
tal data. 

The absence of data at the vicinity of origin (Fig. [6]) 
definitely calls for further experimental investigation of 
diamagnetic instability of strongly correlated electron 
systems. 

One of the properties of the family of curves 6b = 
6b{6x, n) Fig. [S] consists in its universality in a sense that 
the curves being a graphical representation of the rela- 
tionship between the measured properties of IDP (or- 
der parameter 6b and range of existence of IDP 6x) 
and characteristics of the shape of the sample (demag- 
netizing coefficient n) do not depend on concrete Fermi 
surface. Therefore, the data related to different sys- 
tems exhibiting the diamagnetic instability can belong 
to the same curve. One can conclude that plotted in the 
{6b, A) plane where A is the range of existence of IDP 
in terms of internal magnetic field, the family of curves 



(Fig. [6]) gathers onto a single universal, or parent curve 
6b — 6b{6x,n = 0). This is a consequence of the univer- 
sality of the bifurcation theory in description of diamag- 
netic instability of the electron system in a joint space 
(b,a,x) of phase and control variables. Another feature of 
the dependencies represented in Fig.[S]consists in a possi- 
bility of calculation the value of a for given experimental 
arrangement. As it follows from Eq. a straight hori- 
zontal line in Fig. [6] crosses the family of curves in points 
of equal a. Mapping the points belonging to different 
curves onto the parent curve 6b = 6b{A) parameterized 
by a, with careful tabulation provides the values of a 
which can be compared with corresponding values cal- 
culated in different models of Fermi surface. Some ex- 
amples of the mapping of the experimental points onto 
parent curve are shown in Fig. [5] by dash dot lines. It is 
follows from the analysis of the function 6b = 6b{6x, n) 
that predicted dependency on the shape of the sample 
can be detectable for moderate values of reduced ampli- 
tude of dHvA oscillations a w 1.3-3. For greater values 
of a the high precision experiments are needed due to the 
small difference between the curves. 



IV. CONCLUSIONS 

The diamagnetic instability and stable equilibria of the 
system of the strongly correlated electron gas at high 
magnetic field and low temperature are studied by the 
tools of catastrophe theory. While at the center of dHvA 
period the symmetric supercritical pitchfork bifurcation 
gives rise to the second order phase transition on tem- 
perature, the deviation of magnetic field from the value 
corresponding to the center of dHvA period results in 
the change of the phase transition order from the second 
to the weakly first in the nearest vicinity of triple de- 
generate point and first at the rest part of dHvA period 
both on temperature and magnetic field. We calculate 
the generalized condition for DPT occurrence. 

The results of our investigation show a reasonable 
agreement with the experiments. We believe also that 
the observation of DPT and interpretation of data can 
be complicated by the intrinsic structure of IDP which 
can realized through band patterns, bubbles and so on, 
but this is beyond the scope of the present publication. 

We hope that our studies provide consistent framework 
in consideration of diamagnetic instability of strongly 
correlated electron gas, will help in making a quantitative 
systematic analysis of experimental data and stimulate 
the further experiments on investigation of the diamag- 
netic instability. 
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